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We show numerically that any of the constant mean curvature tori first found by 
00 ■ Wente must have index at least eight. 12 
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O ■ 1 Introduction 
Q 

The Hopf conjecture asked if all closed surfaces immersed in M 3 with constant mean curvature 
H must be round spheres. It was proven true when either the surface has genus zero by Hopf 
himself [H], or the immersion is actually an embedding by Alexandrov [H]. However, it does 
not hold in general, and the first counterexamples, of genus 1, were found by Wente [We]. 
Abresch [A] and Walter [Wa] made more explicit descriptions for these surfaces of Wente, 
■ which all have one family of planar curvature lines [Sp]. We call these surfaces the original 
0\ ! Wente tori. 

Constant mean curvature surfaces are critical for area, but not necessarily area mini- 
mizing, for all compactly supported volume-preserving variations. Hence the index - loosely 
\Q ■ speaking, the dimension of area-reducing volume-preserving variations, to be defined in Sec- 
tion 3 - can be positive. If it is zero, the surface is stable. Do Carmo and Peng [CP] showed 
that the only complete stable minimal surface is a plane. Fischer- Colbrie [FC] showed that 
a complete minimal surface in M. 3 has finite index if and only if it has finite total curvature, 
and that the catenoid and Enneper's surface have index 1. Likewise, for surfaces with con- 
stant mean curvature H ^ 0, Barbosa and Do Carmo [BC] showed that only round spheres 
are stable, and Lopez and Ros [LR] and Silveira [Si] independently showed that they have 
finite index if and only if they are compact. This leaves open the question of whether there 
exist surfaces with constant mean curvature H ^ and low positive index, for example with 
index 1. 

The third author [Rl], [R2] showed numerically that the most natural candidates for 
unstable surfaces of constant mean curvature H ^ with low index - the original Wente 
tori - all have index at least 7, and with a numerical experiment suggested that the sharpest 
lower bound is either 8, 9, or 10, and is most likely 9. This leads one to conjecture that all 
closed surface with constant mean curvature H ^ have index at least 9. 

The purpose of this article is to show that the original Wente tori all have index at 
least 8, improving the lower bound of [Rl], [R2]. The final part of our argument relies on 
numerics. 
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2 The original Wente tori 



In this section we shall give a brief description of the original Wente tori, based on [Wa]. 
Later, we shall assume that the mean curvature H is 1/2, but in this and the next section 
we shall only assume that if is a nonzero constant. 

Let ?£ : C/T — > M 3 be a conformal immersion of class C°° where C/T is a compact 
2-dimensional torus determined by the 2-dimensional lattice T. Note that (x, y) are then 
isothermal coordinates on C/T. The fundamental forms and the Gauss and mean curvature 
functions are 

/ = E{dx 2 + dy 2 ), II = Ldx 2 + IMdxdy + Ndy 2 , K = ^=^, H = ^f. 

Since H is constant, the Hopf differential <£> dz 2 is holomorphic, where $ = \{L — N) — %M 
and z = x + iy. Thus $ is constant and ?£ has no umbilics points. Moreover, by a change of 
the coordinates (x, y), we may assume $ = 1 and so M = 0, L = e F + 1, JV = e F - 1, and 
(x, y) become curvature line parameters, where F : C/T — > R is defined by HE = e F . We 
have the equations of Gauss and Weingarten: 

(!) I 1 , _ 

&xy = -^FyX x + -F x ?6 y , M x = H(l + e )?6 X , A/" y = H(l — e )P6 y 



(2) AF + AH sinh F = , 

where A = ^ + and A/" : C/T — > R 3 is the unit normal vector field, i.e. the Gauss 

map. Therefore the problem of finding constant mean curvature immersed tori in M 3 reduces 
to solving the PDE system (1) and (2) by real analytic functions F,J\f, ?6 defined on R 2 and 
doubly periodic with respect to some fundamental lattice rcl 2 . 

In the case of the original Wente tori, in Walter's notation, the solution F of (2) is: 

(3) tanh ( ^ ) = 7 • 7 • cn k (ax) cn^ (ay) , 



where cn& denotes the Jacobi amplitudinus cosinus function with modulus k, and k = 
sin 0, k — sin 0, for 0,0 e (0, vr/2) and + < vr/2, and 

7 = V^T0, 7 = « a = v^H^, a = . 

Lemma 1 ([A], [Wa]). T/ie se^ o/a// original Wente tori are in a one-to-one correspondence 
with the set of reduced fractions £/n G (1, 2). 

For each £/n, we call the corresponding Wente torus VV^/ n . Following Walter's notation, 
each We/ n has either one or two planar geodesic loops in the central symmetry plane: two 
loops if £ is odd, and one loop if £ is even. Each loop can be partitioned into 2n congruent 
curve segments, and £ is the total winding order of the Gauss map along each loop. 

The conditions for double periodicity of the position vector function ?6 are expressed 
in terms of and 0. Walter determined that there is exactly one 

65.354955° 
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that solves one period problem. The other period problem is solved with the correct choice 
of 9 G (0, (tt/2) — 9), and, for any £/n G (1, 2), this correct choice is the unique solution 9 of 



(4) 



71/2 1 + tan 9 tan 9 cos 2 ip dip £ ir / sin 29 



o l-tan^ tan 9 cos 2 </? y/l - sin 2 6» sin 2 y9 n 2 y sin 2(0 + 6>) 



For any £/n G" (1, 2), there is no solution 9 G (0, (7r/2) — 9) of (4). In Table 2 we give some 
values of 9 with respect to £/n. 

Now, if (resp. yg n ) denotes the length of the period of cnk(ax) (resp. cnj(ay)), then 
we have the following lemma: 

Lemma 2 ([Wa]). ?6 : C/r — > Wi/ n C I 3 is a conformal immersion {W^/ n denotes the 
image of ?&), where 

T = span z {(nxi n , 0), (0,yi n )} when£ is odd, and 



T = span z | py^, ^Y^j , (0, w/ien £ is even. 

The curves {[xo,y]\xo = constant} are mapped by ?6 to planar curvature lines ofWi/ n . 
The lengths xe n and y^ n can be computed as follows: 

4 r /2 ^ 4 r /2 

a Jo a/1 - k 2 sin « -A) ^If^ 



3 The definition of index and preliminary results 

The Jacobi operator associated to W</„ is -A/- |//| 2 on C/r, with |//| 2 = £~ 2 (L 2 + 2M 2 + 
N 2 ) = 2H 2 (1 + e~ 2F ) and A/ the Laplace-Beltrami operator associated to the metric I. The 
corresponding quadratic form is 

(6) Q(u,u) = / uC(u)dxdy, 

Jc/t 



where 

Cu = -Au - Vu with V = AH cosh(F) 

and A the Euclidean Laplacian. Note that in equation (6), we are integrating with respect 
to the flat metric on C/r. 

Consider a smooth volume-preserving variation ?6 t of the immersion ?6 with parameter 
t so that is the surface Wi/ n . By reparametrizing the surfaces of the variation, we may 
assume that the variation vector field at t = is uAf for some u G C°°(C/T). Then 

J^area(/¥ t )| t=0 = and J^area^t) ^ = Q(u,u). 
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Furthermore, the volume-preserving condition implies J c ^ r udA = 0. Thus, if 

V = |tj G C°°(C/r); J c/r udA = o}, 
then we can give the following definition (see [BC]): 

Definition 1 We define Ind(<¥(C/r)), the index of the immersion ?£ of C/r, to be the 
maximum of the dimensions of the subspaces of V restricted to which Q is negative definite. 

Since the first derivative of area is zero, and the second derivative is Q(u, u), the index 
in a sense measures the amount of area-reducing volume-preserving variations. 

Let L 2 = L 2 (C/r) = {u e C°°(C/r)| f c/r u 2 dxdy < 00} provided with the inner 

product (u, v)l,2 = f c ^ r uv dx dy. It follows from the standard spectral theorem that the 
operator C = —A — V on C/r has a discrete spectrum of eigenvalues 

A < < ■ ■ ■ / 00 

and corresponding eigenfunctions 

^ 2 ... e c°°(c/r), 

which form an ortho normal basis for L 2 . Moreover, we have the following variational char- 
acterization for the eigenvalues: 

Pi = inf v . (sup 06F . ; mL2=1 f c/r (f)£(l)dxdy}, 

where Vj runs through all j dimensional subspaces of C°°(C/r). 

Lemma 3 ([Rl], [R2]). If C hash negative eigenvalues, then Ind(W / £/ n ) is either k ork — 1. 
Furthermore, if there exists a subspace S C L 2 such that S C C°°(C/r) and dim(<S) = k 
and Q restricted to S is negative definite, then Ind(W^/ n ) > k — 1. 



By Lemma 3, our goal becomes to compute the number of negative eigenvalues of C. 
Now, we use a convenient fact: For the flat torus C/r, with T = span z {(ai, 02), (61, b 2 )}, 
the complete set of eigenvalues of — Aui = onUi are 

( ai & 2 4 -a 2 &i) 2 (( m2&2 ~ m i fl 2) 2 + (miai - m 2 6i) 2 ), 
with corresponding orthonormal eigenfunctions 

c mi ,m 2 ■ (sin or cos) (^r^-((™2&2 - m x a 2 )x + {m x a x - m 2 b 1 )y)^j, 

where m 1 ,m 2 G Z, c mum2 = y/2/(aib 2 - a 2 bi) if \mi \ + \m 2 \ > 0, c , = y/l/(aib 2 - a 2 b 1 ). 
With the aid of Lemma 2 we list 17 of the ccj and U{ in Table 1. 

With the orderings for the eigenvalues as chosen in Table 1, we do not necessarily have 
aii < aij for i < j. However, we still have at /* 00 as i S 00. Choose cx Pi/n (i), ct Pe/n (2), • • • 
the complete set of eigenvalues with multiplicity 1 of the operator —A on the flat torus C/r 
reordered by the permutation py n of N so that a Pj , /n (i) < a p e/n (2) 00. 
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eigenvalues 
oti for t odd 



eigenf unctions 
Ui for t odd 



eigenvalues 
OLi for I even 



eigenfunctions 
Ui for £ even 



ci\ = 



Mi = 



\/nx tn y in 



a\ — 



Ml = 



sin(27ra:/(raz < ; ra )) 
V nx tn y ln /2 



^Jnx ln y ln 



sui(A-Kx/{nxt n )) 
yj nx en y en /4: 



«2 



4tt 2 



«2 



«2 



16tt 2 

fa 



«2 



cos(27ra:/(nz <ra )) 
V nx ln y ln /2 



cos(Airx/(nx en )) 
yj nx£ n yj n /A 



4tT 2 



«3 



aq 



I67H 



«3 



"fa ' % 



a4 



4vr 2 



U4 



sm(2ny/y erl ) 
\/_ nx e „yi n /2 



«4 



An- 

n 2 x'j 



- _|_ 4tT 2 
S'fa 



«4 



sin( 



y/ nx tn y erl /4 



V "fa Wfa ; 



as = 



4tH 



«5 = 



cos(2-ny/y ln ) 
yJnxi n y ln /2 



a§ = 



4tH 



4tt^ 



nx ln y in /< 



V£n ' 



a 6 



167T1 
2 2 



U 6 



sm(Awx/(nx£ n )) 
yj nx tn y ln /2 



a 6 



Att 2 - 

— 9 — T~ 



AlT 2 
^fa 



«6 



sin( 



2 — :r 
"fa 



V nx en y en /4: 



my 
ye 



tt7 = 



16tt 2 



Ui = 



cos(Airx/(nx en )) 



(Xi = 



Att 2 



nx en y in /z 



Att 2 



cos(- 



U7 = 



\/nxi n y ln /A 



a 8 



An' 



. Att 2 

n2x L ' Vtn 



+ 



U 8 



sin(- 

v 7 



yj nxi n yi n /2 



a 8 



v °fa yen ' 

yj nx en y lrl / 



16tt 2 

yen 



u 8 



sm(ATry/y en ) 
yj nxenyen/4: 



a 9 



An 2 



Att 2 



Ug 



a 9 



1671 



Ug 



- « 3 ' 



cos(47n//yi ra ) 
yJnxi n y tn /A 



Att 2 



+ 



Att 2 



cos( 



2tyx 
nx ln 



yj nx tn y ln /2 



"10 



64tt 2 



«10 



my ^ 

ve n ' 



sui(Sttx / (nx ( n )) 
yj nx en y ln /A 



CKn = 



4tt 2 



Att 2 



cos(- 



Mn = 



yJnx en y ln /2 



an = 



64tt< 



Mn = 



cos(87T3:/(ra^ n )) 
yJnx en y ln /A 



yJnx en y ln /A 



Ot\2 



I67H 
2/fa 



«12 



sin(47ry/y£ ra ) 
\/nx tn y trl /2 



"12 



36tt 2 
2 2 



4tt 2 
2/fa 



«12 



c fa ' ^fa ; 



; 13 



I67H 



COs(47T1//l/fa) 

\/ nxt n ye n /2 



'13 



36tt 2 



4tt 2 



cos(- 



«13 



\/ nxi n yt n lA. 



"14 



36?r 2 

" 2 *fa 



Ml4 



sin(67ra;/nx^ n ) 
yJnx tn y ln /2 



"14 



36tt 2 
9 2 
" x fa 



4vr 2 
^fa 



Ml4 



sin(- 



y/nx tn y ln /A 



v n3: <n ^fa . 

it 



"15 = 



36tt 2 



«15 = 



cos(6nx/nx£ n ) 
\/nx tn y ln /2 



am = 



36tt 2 



4vr 2 



«15 = 



nx en ' ^fa ; 
yj nxe n yen/2 



nx in y ln /< 



1 47ry s 



"16 



16tt 2 _|_ 4tt 2 
^fa 



2 2 



«16 



sin( 



"16 



16tt 2 

— 9 — ^~ 

n 2 x 2 n 



16tt 2 
^fa~ 



«16 



sin(- 



iT 

cos( — 



nx ln yen ' 



y/nx en ye n /A 



cos(^2L + ^T 

V "<n !<fa ^ 

yJnx in y ln /A 



"17 



16tt 2 _|_ 4tt 2 



Ml 7 



y/nx en y in /2 



"17 



16tt 2 1 16tt 2 
T 



"-fa 



yfn 



Uyj 



Table 1: The first 17 eigenvalues and eigenfunctions of — A«j = a^. 

The first of the following two lemmas follows from the variational characterization for 
eigenvalues, and the second follows from Lemma 3, the Courant nodal domain theorem, and 
geometric properties of the surfaces Wi/ n : 

Lemma 4 ([Rl], [R2]). Choose u e Z + so that a pe/M < AH. Then Ind(W^ /n ) > u - 1. 

Lemma 5 ([Rl], [R2]). For all n G Z + ; n > 2 we have that Ind(W^/ n ) > In - 2 if t is odd, 
and lnd(Wi/ n ) > n — 2 if £ is even. 

4 The lower bound 8 for Ind(V^/ n ) 



We now show the following: 
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Numerical Result: Ind(W^/ n ) > 8 for all i/n. 

Observe that, although the eigenvalues of C depend on the choice of H, the number of 
negative eigenvalues is independent of H . So without loss of generality we fix H = 1/2. 

By Lemma 5, Ind(W^/ n ) can be less than 8 only if i/n is one of 3/2, 4/3, 5/3, 5/4, 7/4, 
6/5, 8/5, 8/7, 10/7, 12/7, 10/9, 14/9, or 16/9. Lemma 4 also gives explicit lower bounds for 
the index, since we know the values of xe n and yg n numerically by formula (5), and hence 
we know the ci Pe/n (i) (see Table 1). Lemma 4 implies that the index is at least 8 when i/n is 
5/4, 6/5, 8/7, 10/7, or 10/9. Thus we only need to consider the following eight surfaces: 

W3/2, W4/3, W5/3, W7/4, W 8 /5, W12/7, Wu/g, and W 16/9 . 

For these surfaces we list, in Table 2, the corresponding 9, xe n , y^ n and lower bounds 
for index. These approximate values for 9, xe n , and y^ n were computed numerically using 
formulas (4) and (5) and the software Mathematica. Recall that always 9 = 65.354°. 



w t/n 


9 




Ven 


Lemma 4 
lower bound 
for Ind(Wi n ) 


Lemma 5 
lower bound 
for lnd(W tn ) 


W 3/2 


17.7324° 


2.5556 


4.2131 


2 


2 


w 4/3 


12.7898° 


3.2767 


6.3355 


6 


1 


w 5/3 


21.4807° 


1.7557 


2.6402 


2 


4 


w 7/4 


22.8449° 


1.3315 


1.9447 


2 


6 


w 8/5 


20.1374° 


2.0842 


3.2321 


2 


3 


w 12/7 


22.3044° 


1.5150 


2.2380 


2 


5 


w l4/9 


19.1243° 


2.2970 


3.6514 


4 


7 


w l6/9 


23.2182° 


1.1872 


1.7208 


2 
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Table 2: xe n , yen are computed using the value H = 1/2. 

We will find specific spaces on which C is negative definite, for these eight surfaces. 

Let N be an arbitrary positive integer. Consider a finite subset {Hi — u^, . . . ,un = Ui N } 
of the eigenfunctions Ui of —A on C/r, defined in Section 3, with corresponding eigenvalues 
a.j — aj., j = 1, . . . , N. If we consider any u = J2iLi a iUi e span{iti, . . . ,%}, a±, . . . , £ 
R, then J c ^ r uC(u)dx dy = Y^Hj=i a i{^j^i3 ~ hj) a j-, where 6^ := J c ^ r VuiUj dxdy. So we 
have J c ^ r uC(u)dx dy < for all nonzero u G span{-u 1? . . . , u N } if and only if the matrix 
(aijSij — fcjj)jj=i,...,jv is negative definite. Lemma 3 then implies: 

Theorem 1 ([Rl], [R2]). If the N x N matrix (aijSij — j=i at is negative definite, then 
lnd(W e/n )>N-l. 



We, n 


Ml 


u 2 


U3 


M 4 


M 5 


M 6 


M 7 


u 8 


M 9 


W3/2 


Ml 


u 2 


u 3 


M 4 


M 5 


M 7 


M 8 


M 9 


M17 


^4/3 


Ml 


u 2 


U3 


M4 


M 5 


M 6 


M 7 


M 8 


M 9 


^5/3 


Mi 


u 2 


U3 


M 5 


M 6 


My 


M 8 


M 9 


U15 


W7/4 


Mi 


u 2 


U3 


M 6 


M 7 


M 8 


M 9 


M14 


U15 


^8/5 


Mi 


u 2 


U3 


M4 


M 5 


«10 


Mn 


U12 


Uis 


Wl2/7 


Mi 


u 2 


U3 


M4 


M 5 


«10 


Mn 


U12 


Uis 


W^i4/9 


Mi 


u 2 


U3 


M4 


M 5 


u w 


Mn 


Ul2 


Uis 


Wl6/9 


Mi 


u 2 


Us 


M 4 


M 5 


UlO 


Mn 


U12 


U13 
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Table 3: Eigenj 'unctions of —A producing 9 -dimensional spaces on which Q is negative 

definite. 

Definition 2 Given A, B even integers and £,n G Z + , we now define the following basic 
integrals: 

W,n,A,B) = -±- f ^J^v(cos^) A (cos^Ydvdx, 

h <*. n) = ^ sr Jl r /4 v co S (s) *«**■ 
w r " i/4 r /4 " «» (s) d **> 

W ") = /r" ,/4 /o""' 4 " cos (tSw) d vdx, 

n) = — §— C XtJ4 r n/4 V cos f ^ cos ( dydx, 

*V 5 / nx tn y en JO JO \nx en J \yt„ J w ' 

J 5 (£, n ) = —5— C XeJi ( n W4 V cos ( ^) cos f dydx, 

° v ' y nx ln y ln JO JO \nx ln ) \nx in J & ' 

w. ") = =£= r'- /4 jt /4 17 «» (S) ™ (e) <W». 

/ 7 (C, „) = f ■"*>•>" [«'■" V cos ( i^) cos ( in) dydx. 

'V ' ' m(n!»n J J \nx in J \ye. n J a 

Now, for each surface W£/ n given in Table 2, we will fix N = 9 and choose the subset 
{Hi, . . . , -ug} such that the matrix (atjSij — &ij)ij=i,...,7v is negative definite. These choices are 
given in Table 3. With these choices for Ui, we have the following lemma: 

Lemma 6 With the choices given in Table 3, all elements of the eight matrices M.{£,n) := 
(aijSij — &ij)ij=i,...,9 can be expressed in terms of the basic integrals Io(£, n, A, B) and Ij(£, n) 
for A, B even and j = 1, 2, . . . , 7. 

Proof: The symmetries V(x, y) = V(—x, y) = V(x, —y) = V {^f- — x, y) = V (x, ^ — y) 
of V and the identities cos(a ± b) — cos(a) cos(6) =F sin(a) sin(fe), sin(a ± b) — sin(a) cos(6) ± 
sin(6) cos(a) give the relations shown in Table 4, proving the lemma. 
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For 



W, 



3/2 



Mi,i = -32 J (3, 2, 0, 0), M 4>4 = a 4 - 64(/ (3, 2, 0, 0) - / (3, 2, 0, 2)) 

M 5i5 = a 5 - 64 J (3, 2, 0, 2), M 7>7 = a 7 - 64 J (3, 2, 2, 0) 

M 17 ,n = «i 7 - 64(J (3, 2, 0, 0) - J (3, 2, 0, 2) - J (3, 2, 2, 0) + 2 J (3, 2, 2, 2)) 

M M = ai - 32/ (3, 2, 0, 0), for i = 2, 3, 8, 9 

Mjj = for i < j with i,j £ {1,2,3,4,5,7,8,9, 17}. 



For 

W 4/3 



Mi,i = -48 J (4, 3, 0, 0), M 2i2 = a 2 - 48 J (4, 3, 0, 0) + 2 J 2 (4, 3) 

M 3i3 = a 3 - 48 J (4, 3, 0, 0) - 2 J 2 (4, 3), M 4 , 4 = a 4 - 48 J (4, 3, 0, 0) + 2 J 4 (4, 3) 

M 5)5 = a 5 - 48 J (4, 3, 0, 0) - 2 J 4 (4, 3), M 6>6 = a 6 - 48 J (4, 3, 0, 0) + 2 J 4 (4, 3) 

M 7j7 = a 7 - 48 J (4, 3, 0, 0) - 2 / 4 (4, 3) 

M 8i8 = a 8 - 384 (J (4, 3, 0, 2) - J (4, 3, 0, 4)) 

M 9i9 = a 9 - 96 (4 J (4, 3, 0, 4) - 4 J (4, 3, 0, 2) + / (4, 3, 0, 0)) 

Mi, 3 = -2v^/i(4, 3), Mi, g = -48v^(-/ (4, 3, 0, 0) + 2/ (4, 3, 0, 2)) 

M 4>6 = -48(-/ (4, 3, 0, 0) + 2/o(4, 3, 0, 2)) + 2^(4, 3) 

M 5 , 7 = -48(-/ (4, 3, 0, 0) + 2/o(4, 3, 0, 2)) - 2/^4, 3) 

M 3)9 = —4 / 4 (4, 3), all other M i:j with % < j < 9 are zero. 



For 



W, 



5/3 



Mi,i = -48 / (5, 3, 0, 0), M 2>2 = a 2 - 48 / (5, 3, 0, 0) + 2 A (5, 3) 

M 3i3 = « 3 - 48 / (5, 3, 0, 0) - 2 A (5, 3), M 5 , 5 = a 5 - 96 / (5, 3, 0, 2) 

M 6)6 = a 6 - 48 / (5, 3, 0, 0) + 2 / 2 (5, 3), M 7J = a 7 - 48 / (5, 3, 0, 0) - 2 / 2 (5, 3) 

M 8i8 = a 8 - 48 / (5, 3, 0, 0) + 2 / 4 (5, 3), M 9 , 9 = a 9 - 48 / (5, 3, 0, 0) - 2 / 4 (5, 3) 

Mi5,i5 = «i5 - 96 / (5, 3, 2, 0), M 1)7 = -2^(5, 3) 

-^3,15 = — 2(/i(5, 3) + / 2 (5, 3)), all other pertinent with 2 < j are zero. 



For 

W7/4 



M M = ai - 64 / (7, 4, 0, 0) for % = 1, 2, 3, 6, 7, 8, 9, 14, 15 
all other pertinent Mj ~ with i < j are zero. 



For 

W 8 /5, 
W12/7, 
W14/9, 
Wl6/9 



M^i = -16n J (4 n, 0, 0), M 2>2 = a 2 - 16n / (£, n, 0, 0) + 2h(£, n) 

M 3 [ 3 = a 3 -16nl (£,n,0,0)-21 1 (£,n), M 4>4 = a 4 -16n I (£, n, 0, 0) + 2/ 4 (£, n) 

M 5j5 = a 5 - 16n J (4 ra, 0, 0) - 2/ 4 (£, n) 

M 10 ,i = aio - 16n J (4 n, 0, 0) + 2/ 3 (£, n) 

M 11)U = a n - 16n / (£, n, 0, 0) - 2I 3 (£, n) 

M l2 ,i2 = ai2 - 16n J (4 n, 0, 0) + 2I 7 (£, n) 

M 13) i 3 = ai 3 - 16n / (£, n, 0, 0) - 2I 7 (£, n), M 1)3 = -2yfih{i, n) 
M 1>u = -2V2I 2 (£,n), M 2)10 = -4 J^, n) + 4 / 5 (£, n) 
M 3 , n = -4 / 5 (£,n), M 4 , 12 = -2h(£,n) + 2I 6 (£,n) 

-^5,13 — —^h(£,n) — 2I 6 (£,n), all other pertinent Mjj with i < j are zero. 



Table 4: Elements M it j of the symmetric matrices A4(£,n) expressed in terms of the 
basic integrals. We have chosen here to index the using the counters associated to otj 

and Uj, rather than ctj and Uj. 

By numerical methods, we can estimate that all of the relevant Ij (£, n) for j > 1 are 
approximately zero, and that 

/ (3, 2, 0, 0) 0.2968, / (3, 2, 2, 0) 0.2304, / (3, 2, 0, 2) = 0.2408, / (3, 2, 2, 2) = 0.1947, 
/ (4,3,0,0) 0.1077, / (4,3,0,2) = 0.0776, / (4,3,0,4) = 0.0667, / (5,3,0,0) 0.4532, 
/ (5,3,2,0) 0.3910, / (5,3,0,2) = 0.4046, / (7,4,0,0) 0.6072, / (8, 5, 0, 0) 0.1878, 
/ (12, 7, 0, 0) 0.2652, / (14, 9, 0, 0) 0.0841, / (16, 9, 0, 0) 0.3419. 
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These values were computed with a Mathematica program using the NIntegrate and Jaco- 
biCN commands, and the program is available at the web site of the third author. One note 
of warning is that Mathematica has different conventions than Walter's paper, and hence 
crik in [Wa] is equivalent to cn^ in Mathematica. We include a sample of our code in the 
Appendix. 

Now we can make approximations for the eight matrices Ai(£,n). 
The matrix .M(3, 2) is approximately 



M(S,2) 



( - 



9.50 





























-7.99 





























-7.99 





























-1.36 





























-13.2 





























-8.70 





























-5.76 





























-5.76 





























-5.50 



\ 



and all nondiagonal terms are known to be zero by rigorous mathematical computation, and 
all nonzero entries have been computed only numerically. 
.M(4,3) is approximately the nondiagonal matrix 



/ 



jM(4,3) 



5.17 





O 

















-3.23 





-3.53 























O 





-3.53 

















o 











-3.78 





-2.29 























-3.78 





-2.29 

















-2.29 





-3.78 























-2.29 





-3.78 





























-0.25 





3.23 





o 

















-2.21 



and again here all entries that are have been computed mathematically rigorously, and all 
nonzero entries have been computed only numerically. The symbol O denotes an entry that 
has been computed numerically to be approximately zero, but not mathematically rigorously. 
We shall continue to use these conventions in all remaining matrices. 
.M(5,3) is approximately the diagonal matrix 



M(5,3) 



21.8 














o 














-20.3 





























-20.3 

















o 











-33.2 





























-16.1 





























-16.1 





























-14.7 





























-14.7 











o 

















-24.7 



\ 
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.M(7, 4) is approximately the diagonal matrix 



-M(7,4) 



38.9 





























-37.5 





























-37.5 





























-33.3 





























-33.3 





























-27.0 





























-27.0 





























-26.3 





























-2e 



Ai(8,5) is approximately the diagonal matrix 



A<(8,5)« 



/ 


-15.0 





o 











o 





\ 





-13.6 











o 













o 





-13.6 











o 



















-10.9 











o 



















-10.9 











o 







o 











-9.2 













o 





o 











-9.2 



















o 











-8.0 





V 


























-8.0 / 



(12, 7) is approximately the diagonal matrix 



A4(12,7)« 



/ 


-29.7 

















O 





\ 







-28.3 











o 













O 





-28.3 











O 



















-21.5 











o 



















-21.5 











o 







O 











-24.1 













o 





O 











-24.1 



















O 











-18.7 





V 














o 











-18.7 / 



.M(14, 9) is approximately the diagonal matrix 



M(U,9) 



( 


-12.1 





o 



























-11.7 











o 















o 





-11.7 











o 





















-9.1 











o 





















-9.1 











o 









o 











-10.6 















o 





o 











-10.6 





















o 











-8.3 







V 














o 











-8.3 


/ 
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.M(16, 9) is approximately the diagonal matrix 



M(16, 9) 



/ 


AQ 9 


n 
u 




n 
u 


n 
u 


n en 

U <J 


n 
u 


n 
u 


\ 







-47.9 











O 


















-47.9 








o 





















-35.6 








o 





















-35.6 








o 









O 











-43.7 


















O 








-43.7 





















O 








-32.8 







V 














O 








-32i 





All eight of these matrices are 9x9 and negative definite. Hence Theorem 1 implies 
the numerical result. 

5 Appendix: the Mathematica code 

The following is a Mathematica code for computing the values 7o(4, 3, 0,0), 7o(4, 3, 0,2), 
7 (4,3,0,4), 7i(4,3), J 2 (4,3), J 4 (4,3), and the elements of the matrix A4(4,3). The seven 
other needed codes for different I and n were written similarly. 

H = 1/2; kl = Sin[thetal]; k2 = Sin [theta2] ; 

gammal = Sqrt [Tan [thetal] ] ; gamma2 = Sqrt [Tan [theta2] ] ; 

alphal = Sqrt [4 H Sin [2theta2] /Sin [2 (thetal + theta2)]]; 
alpha2 = Sqrt [4 H Sin [2thetal] /Sin [2 (thetal + theta2)]]; 

F = 4ArcTanh [gammal gamma2 JacobiCN [alphal x, kl~2] JacobiCN [alpha2 y, k2~2]]; 
V = 4 H Cosh[F] ; 

ell =4; n = 3; 

thetal = 2 Pi (12.7898/360); theta2 = 2 Pi (65.354955354/360); 

xO = 3.2767; yO = 6.3355; 

Print ["I_0(4,3,0,0) is " , I0x4c3c0c0x = (l/(n xO y0)) NIntegrate [ 

V , {x, 0, x0/4>, {y, 0, yO/4}]] ; 

Print ["I_0(4,3,0,2) is " , I0x4c3c0c2x = (l/(n xO y0)) NIntegrate [ 

V (Cos [2 Pi x/xO])-0(Cos[2 Pi y/y0])"2,{x, 0, xO/4}, {y, 0, yO/4}]] ; 

Print ["I_0(4,3,0,4) is " , I0x4c3c0c4x = (l/(n xO y0)) NIntegrate [ 

V (Cos [2 Pi x/xO])~0(Cos[2 Pi y/y0])-4,{x, 0, xO/4}, {y, 0, yO/4}]] ; 

Print ["I_l(4,3) is \Ilx4c3x = (8/(n xO y0)) (NIntegrate [ 

V (Cos [4 Pi x/(n xO)]) ,{x, 0, xO/4}, {y, 0, yO/4}] + NIntegrate [ 

V (Cos [4 Pi x/(n x0)]),{x, xO/4, 2 xO/4}, {y, 0, yO/4}] + NIntegrate [ 

V (Cos[4 Pi x/(n x0)]),{x, 2 xO/4, n xO/4}, {y, 0, yO/4}])]; 

Print ["I_2(4,3) is ",I2x4c3x = (8/(n xO y0)) (NIntegrate [ 

V (Cos [8 Pi x/(n xO)] ) ,{x, 0, xO/4}, {y, 0, yO/4}] + NIntegrate [ 

V (Cos [8 Pi x/(n x0)]),{x, xO/4, 2 xO/4}, {y, 0, yO/4}] + NIntegrate [ 

V (Cos[8 Pi x/(n x0)]),{x, 2 xO/4, n xO/4}, {y, 0, yO/4}])]; 
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Print["I_4(4,3) is ",I4x4c3x = (8/(n xO yO)) (NIntegrate [ 

V (Cos [4 Pi x/(n xO)]) (Cos [4 Pi y/yO] ) , {x,0,x0/4} , {y ,0,y0/4}] +NIntegrate [ 

V (Cos [4 Pi x/(n xO)]) (Cos [4 Pi y/yO] ) ,{x,x0/4,2 x0/4> , {y , ,y0/4>] +NIntegrate [ 

V (Cos [4 Pi x/(n xO)]) (Cos [4 Pi y/yO]),{x, 2 xO/4, n xO/4}, {y, 0, yO/4}])]; 



aa = 0; bb = 0; 



alphal = aa(4 N [Pi"2] / (n~2 x0~2)) + bb (4 N [Pi~2] / (y0~2) ) ; 



aa =4; bb = 0; alpha2 = aa(4 N [Pi~2] / (n~2 x(T2)) + bb (4 N [Pi~2] / (y0~2) ) ; 

alpha3 = alpha2; 



aa = 1; bb = 1; alpha4 
alpha5 = alpha4; alpha6 



aa(4 N[Pi~2]/(n~2 x(T2)) + bb (4 N [Pi~2] / (y0~2) ) ; 
alpha4; alpha7 = alpha4; 



aa = 0; bb = 4; 
alpha9 = alpha8; 



alpha8 = aa(4 N [Pi~2] / (n~2 x(T2)) + bb (4 N [Pi~2] / (y0~2) ) ; 



Print 
Print 
Print 
Print 
Print 
Print 
Print 
Print 
Print 
Print 
Print 
Print 



["M(l,l) 
["M(2,2) 
["M(3,3) 
["M(4,4) 
["M(5,5) 
["M(6,6) 
["M(7,7) 
["M(8,8) 
["M(9,9) 
["M(l,9) 
["M(4,6) 
["M(5,7) 



is 
is 
is 
is 
is 
is 
is 
is 
is 
is 
is 
is 



-48 I0x4c3c0c0x] ; 
alpha2 - 48 I0x4c3c0c0x] 
alpha3 - 48 I0x4c3c0c0x] 
alpha4 - 48 I0x4c3c0c0x] 
alpha5 - 48 I0x4c3c0c0x] 
alpha6 - 48 I0x4c3c0c0x] 
alpha7 - 48 I0x4c3c0c0x] 
alpha8 - 384 (I0x4c3c0c2x - I0x4c3c0c4x)] ; 
alpha9 - 96 (I0x4c3c0c0x + 4 I0x4c3c0c4x - 4 I0x4c3c0c2x)] 
-48 N[Sqrt[2]] (-I0x4c3c0c0x + 2 I0x4c3c0c2x)] ; 
-48 (-I0x4c3c0c0x + 2 I0x4c3c0c2x)] ; 
-48 (-I0x4c3c0c0x + 2 I0x4c3c0c2x)] ; 
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